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Abstract 

A Lyapunov-based approach for the trajectory generation of an A^-dimensional Schrodinger 
equation in whole M. N is proposed. For the case of a quantum particle in an A^-dimensional 
decaying potential the convergence is precisely analyzed. The free system admitting a 
mixed spectrum, the dispersion through the absolutely continuous part is the main obstacle 
to ensure such a stabilization result. Whenever, the system is completely initialized in the 
discrete part of the spectrum, a Lyapunov strategy encoding both the distance with respect 
to the target state and the penalization of the passage through the continuous part of the 
spectrum, ensures the approximate stabilization. 



1 Introduction 



1.1 Main results 



We consider a quantum particle in an iV-dimensional space, with a potential V(x), 
and coupled to an external (laser) field 1 1— > u(t) & R through its dipole moment 
ju(x). Under appropriate change of scales, the system's wavefunction evolves fol- 
lowing the Schrodinger equation 

i^-(t,x) = -A x ¥(t,x) + (V(x)-u(t)ii(x)) i ¥(t,x) : x G R. , (1) 
at 

W(0,x)=W (x). (2) 
This is a bilinear control system, denoted by (E), where 
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• the control is the external field u : M+ — > M, 

• the state is the wave function *P : R + x -> C with *P(f) G § for every £ > 0, 
with§ = {<pGL 2 (M A ';C) | ||<p|| L 2 = 1}. 

We distinguish between four different situations: dimension N > 4; dimension N = 
3; dimension N = 2; and dimension N = I. For each of these cases, we will assume 
some appropriate decay assumptions for the potential V(x). Indeed, through this 
paper, we will assume the following assumption: 

Decay assumption (A) We assume for the potential V that zero is neither an eigen- 
value nor a resonance of the Hamiltonian Hq = — A + V. Furthermore, we as- 
sume one of the following assumptions (depending on the space dimension N) 

• N=l: (l + \x\)Vel}(R) m; 

• N = 2: \V(x)\<C(l + \x\)- 3 - £ EH; 

• N = 3: V G d- £ (R 3 ) nzJ+ e (M 3 ) EH; 

• N > 4: V G L 1 and (1 + |x| 2 ) 7 / 2 V(x) is abounded operator on the sobolev space 
H v for some v > and y > n + 4 ||29l . 

A brief discussion on the origin of the above assumption on the potential V is 
provided in the Subsection ll.3l As one will see these decay assumptions are chosen 
to assure relevant dispersive estimates. 

Furthermore, note that, under the decay assumption (A) on the potential V, the free 
Hamiltonian Hq = — A + V(x) admits a mixed spectrum: 

0(H ) = C7disc(#o) n Oae(H ), 

where the discrete spectrum Odisc(#o) contains a finite number of eigenvalues of 
finite multiplicities and the essential spectrum is actually an absolutely continuous 
spectrum o ac (Ho) = [0,°°). Under the decay assumption (A), this decomposition 
of the spectrum for the ID case is a classical result of the earliest days of quantum 
mechanics (in fact one only needs V G see e.g. 11391 , Sec. XIII.4). For the 

2D case, one can find a proof in ||4~3l . The 3D case has been proven in [23]. Finally, 
the decomposition for the A^-dimensional case, with Af > 3, is a classical result as 
the potential is a short range potential in the sense of Agmon Q]|. 

Concerning the bound states, {0j}^l o , we know that G H 2 (R N ,C). Moreover, 
the decay assumption (A) on the potential implies that V G Lj and V G M\ oc the 
local Stummel class (see @, page 8, for a definition). This ensures the exponential 
decay of the eigenfunctions {0j}jl o (see e.g. [2], page 55, Corollary 4.2). 

Let us recall the following classical existence and uniqueness result for the open- 
loop system (OQ)-©. A proof of this result is given in the Appendix. 

Proposition 1 Let the potential V(x) satisfy the decay assumption (A) and consider 
[X G L°°(R N , R). Let ^ G S, T > and u G C°( [0, T) , R). There exists a unique weak 
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solution of ©-©, i.e. a function Y G C°([0,7],§) nC 1 ([0,r],fl"- 2 (R iV ,C)) smc/? 
that 

+ i I e~ iHo(f -^u(s)n(x) } ¥(5)d5 in L 2 (M. N , C) for t G [0, 71 , (3) 
Jo 

and then © noWs m /7- 2 (M w , C). 

If moreover, G H 2 (M. N , C) an J multiplication by ju(x) defines a bounded oper- 
ator over H 2 (R N ,R), then W is a strong solution, i.e. W G C°([0, T],H 2 (R N ,C)) 
PC 1 ([0, 7],L 2 (M W , C)), ffe; equation © no/ds m L 2 (M. N , C) /or f G [0, 7] and ffo» 
initial condition © no/Js in H 2 (M. N , C). 

7/ze weafc f res/?, strong ) solution is continuous with respect to the initial condition 
for the C°([0, 7],L 2 ) -topology (resp. C°([0, 7],/7 2 ) -topology). 

Assuming the potential V(x) such that the discrete spectrum <7disc(#o) is non- 
empty, we are interested here in stabilizing one of the eigenfunctions in this discrete 
part. Fixing £ > to be a small positive constant and considering <j) to be a normal- 
ized eigenfunction in this discrete part, we are interested in designing a feedback 
law « e ( x I / ) such that, the solution *P(?,x) of ©-© satisfies 

liminf | 0F(f ,x) | <j)(x)) \ 2 > 1 - £. (4) 

t — >oo 

Here 

\ C) = [ MxK(x)dx, 

denotes the Hermitian product of L 2 (M. N ,C). Note that, *J? and <j) living on the unit 
sphere § of L 2 (R N ,C), the limit © denotes the approximate stabilization of the 
eigenfunction § (x) . 

Note that, even though the feedback stabilization of a quantum system necessitates 
more complicated models taking into account the measurement backaction on the 
system (see e.g. H26I25I34I0 . the kind of strategy considered in this paper can be 
helpful for the open-loop control of closed quantum systems. Indeed, one can apply 
the stabilization techniques for the Schrodinger equation in simulation and retrieve 
the control signal that will be then applied in open-loop on the real physical system. 
As it will be detailed below, in the bibliographic overview, such kind of strategy has 
been widely used in the context of finite dimensional quantum systems. 

The main result of this article is the following one. 

Theorem 2 Consider the Schrodinger equation ©- ©. We suppose the potential 
V(x) to satisfy the decay assumption (A) and we take \l G ) nL°°(R ). We as- 

sume the discrete spectrum Odi sc ofH^ = —A + V(x) to be non-empty. We consider 
moreover the following assumptions: 
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Al = If =0 arfj where -}f = 

are different normalized eigenfunctions in the 

discrete spectrum of Ho. 

A2 the coefficient OCq corresponding to the population of the eigenfunction (f>o in 
the initial condition is non-zero: CCq 7^ 0. 

A3 the Hamiltonian Hq admits non-degenerate transitions: Xj x — 7^ Ay 2 — Xk 2 
for (ji,k\) 7^ (72,^2) andwhere {A 7 }jl are different eigenvalues of the Hamil- 
tonian Hq; 

A4 the interaction Hamiltonian /i(x) ensures simple transitions between all eigen- 
functions of Ho: 

(ju^-I^/O Vj7^e{0,l,...,M}. 

Then for any e > 0, there exists a feedback law u{t) = u £ (*H{t)) (that we will con- 
struct explicitly), such that the closed-loop system admits a unique weak solution 
in C°([0, r],S) nC^fO, T],H- 2 (R N ,C))- Moreover the state of the system ends up 
reaching a population more than (1 — e) in the eigenfunction 0o (approximate sta- 
bilization): 

liminf I | fo(x)) | 2 > 1 - £. 

f — >oo 

If moreover multiplication by fl(x) defines a bounded operator over H 2 (R N ), then 
¥ is a strong solution, i.e. ¥ E C°([0,T],H 2 (R* C)) nC 1 ([0,7 , ],L 2 (R iV ,C)). 

Remark 3 In this Theorem Jzf denotes \J p>2 L p (R N ). 

Remark 4 Note that, as the initial state is a linear combination of the bound states, 
we have in particular E H 2 and decays exponentially. 

Remark 5 Note that, here a finite dimensional approximation of the system by re- 
moving the continuous part of the spectrum is not sufficient to treat the stabilization 
problem. In fact, even if the system is initialized in the discrete part of the spectrum 
(as assumed in Al ), the interaction Hamiltonian /1 will make the solution leave this 
discrete part. The state of the system will therefore leave this subspace just after 
the initial time. 

The assumptions Al through A4 can be relaxed significantly. However, as the final 
result with the relaxed assumptions may seem too complicated, we will discuss this 
relaxations, separately, in Section [5l 

1.2 A brief bibliography 

The controllability of a finite dimensional quantum system, i^V = (Ho + u(t) Hi)^ 
where *F E C N and Hq and Hi are N x N Hermitian matrices with coefficients in 
C has been completely explored H46I37I3I4I49I . However, this does not guarantee 
the simplicity of the trajectory generation. Very often the chemists formulate the 
task of the open-loop control as a cost functional to be minimized. Optimal control 
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techniques (see e.g., 114210 and iterative stochastic techniques (e.g, genetic algo- 
rithms 11321 ) are then two classes of approaches which are most commonly used for 
this task. 

When some non-degeneracy assumptions concerning the linearized system are sat- 
isfied, ll35l provides another method based on Lyapunov techniques for generating 
trajectories. The relevance of such a method for the control of chemical models 
has been studied in [361. Since measurement and feedback in quantum systems 
lead to much more complicated models and dynamics than the Schrodinger equa- 
tion [12613411 , the stabilization techniques presented in 11331 are only used for gen- 
erating open-loop control laws. Simulating the closed-loop system, we obtain a 
control signal which can be used in open-loop for the physical system. Such kind 
of strategy has already been applied widely in this framework HI 31451 . 

The situation is much more difficult when we consider an infinite dimensional 
configuration. Concerning the controllability problem, very few results are avail- 
able [148181 1011 . In H8I101 the controllability of a particle in a moving one dimen- 
sional quantum box has been studied. A local controllability result is therefore 
provided using the return method [fT6ll . In |[T2l . applying some geometric control 
tools, the authors provide a quite general result concerning the controllability of 
discrete-spectrum Schrodinger equation. Finally, in 11471 . the authors consider the 
controllability of some particular Schrodinger equations with continuous spectra. 

Concerning the trajectory generation problem for infinite dimensional systems still 
much less results are available. The very few existing literature is mostly based on 
the use of the optimal control techniques H6I71 . The simplicity of the feedback law 
found by the Lyapunov techniques in [1351911 suggests the use of the same approach 
for infinite dimensional configurations. However, an extension of the convergence 
analysis to the PDE configuration is not at all a trivial problem. Indeed, it requires 
the pre-compactness of the closed-loop trajectories, a property that is difficult to 
prove in infinite dimension. This strategy is used, for example in [fT5l . 

Let us mention some strategies for proving the stabilization of infinite dimensional 
control systems. One can try to build a feedback law for which one has a strict 
Lyapunov function. This strategy is used, for example, for hyperbolic systems of 
conservation laws in lfT4ll . for the 2-D incompressible Euler equation in a simply 
connected domain in ifTTl . see also [fT9ll for the multi-connected case. For systems 
having a non controllable linearized system around the equilibrium considered, the 
return method often provides good results, see for example | T6l for controllable 
systems without drift and [12010 for Camassa-Holm equation. In the end, we refer 
to |[T8l for a pedagogical presentation of strategies for the proof of stabilization of 
PDE control systems. 

In this paper, we propose a Lyapunov-based method to approximately stabilize a 
particle in an N dimensional decaying potential under some relevant assumptions. 
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We assume that the system is initialized in the finite dimensional discrete part of the 
spectrum. Then, the idea consists in proposing a Lyapunov function which encodes 
both the distance with respect to the target state and the necessity of remaining 
in the discrete part of the spectrum. In this way, we prevent the possibility of the 
"mass lost phenomenon" at infinity. Finally, applying some dispersive estimates of 
Strichartz type, we ensure the approximate stabilization of an arbitrary eigenfunc- 
tion in the discrete part of the spectrum. 

The ideas of this paper (a short and simplified version is already published as a 
communication Il33l0 have been recently adapted to the case of a quantum particle 
in an infinite potential well IfTTTl . In IfTTI . as we are dealing with a pure discrete 
spectrum, much less restrictive assumptions are needed to ensure the approximate 
stabilization of the system. 

As it can be remarked through the bibliography, except for a very few results [|47ll , 
all the previous work on the control of the infinite dimensional quantum systems 
deal with discrete- spectrum Schrodinger equations. It seems that the techniques of 
this paper and the possibility of the relaxations, explained in Section[5l can open a 
new gateway to investigate this class of quantum systems. 



1.3 Free dynamics and dispersive estimates 



Before treating the control problem, let us have a look at the behavior of the 
system in the absence of the control field (u(t) = 0). We will denote by S(t) = 
exp(-itHo) the Co-semigroup on L (R , C) spanned by the infinitesimal generator 
(— A + V(x))/i. Note in particular that, S(t) induces an isometry overL (R , C): 

\\s(t)w\\ L 2=\\w\y- 

Moreover, we denote by Pdi sc the projection operator over the discrete subspace 
generated by the bound states and defined on L (R , C). Finally, P ac denotes the 
projection over the orthogonal subspace: P ac = Id — Pdisc- 

The discrete part of the freely evolving solution P disc 5(?) l I / o represents a quasi- 
periodic behavior: 

M M 
*P0,disc = Pdisc^O = £ 0Cj$j(x) PdiscSW^O = £ CCje-^Qjix). 

The continuous part, however, represents a dispersive behavior. In this subsection, 
we provide a very brief overview of the dispersive estimates and in particular the 
ones we use in this paper. 
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In heuristic terms, for the potential-free problem V = 0, the explicit solution 



(e ltA W ) (*) = C N t~ N > 2 f ^ ¥ (y)dy, 

implies the dispersive estimate ll44ll 

sup\t\ N / 2 \\e ItA ¥ \\ L ~ {RN) < \\yr\\ Llm V ¥ eL\R N )nL 2 (R N ). 

t>0 

For general V ^ 0, no explicit solutions are available and therefore one needs to 
proceed differently. Consider the perturbed Hamiltonian Hq = — A + V, we seek 
to prove similar estimates on the time evolution operator S(f)P ac = e~' tH °F ac . The 
projection onto the absolutely continuous spectrum of Hq is needed to eliminate 
bound states which do not decay over any length of time. We, therefore, have the 
following dispersive estimate: 

Theorem 6 Under the decay assumption (A) on the potential V, we have 

||S(?)Pa C ||l^oo<M-i (5) 

Such dispersive estimates have a long history. For exponentially decaying poten- 
tials, Rauch [38] proved dispersive bounds in exponentially weighted L 2 -spaces. 
Jensen and Kato ll27l replaced exponential with polynomial decay and obtained 
asymptotic expansions of e~' tH ° (in terms of powers of in the usual weighted 
L 2 '° spaces. The first authors to address a dispersive estimate of the form © were 
Journee, Soffer, and Sogge ll29ll . They were able to prove the dispersive estimate © 
under the fourth case of the decay assumption (A) for N >3. 

Concerning the case N = 3, following a large amount of results H51I40I22H . finally 
Goldberg ETI proved the dispersive estimate © under the third case of the decay 
assumption (A). In contrast, trying to adapt these results to higher dimensions has 
lead Goldberg and Visan E4l to show that for N > 4, © fails unless V has some 
amount of regularity, i.e., decay alone is insufficient for © to hold if Af > 4. 

The one-dimensional case was open until recently. Weder ||50l proved a version 
of Theorem [6] under the stronger assumption that f" M |V(x) | (1 + \x\) 3 ^ 2+£ dx < °°. 
Finally, in a similar way to ll50ll . Goldberg and Schlag Il22l were able to prove © 
under the first case of the decay assumption (A). 

Finally, concerning the two-dimensional case, Yajima ||52l and Jensen, Yajima |[28l 
proved the L P (M?) bounded-ness of the wave operators under stronger decay as- 
sumptions on V(x) (than the decay assumption (A)), but only for 1 < p < °°. Hence 
their result does not imply ©, but LP — > L p estimates for 1 < p < 2. The first paper 
to provide an L 1 — > L°° dispersive estimate of the form © in two dimensions was 
that of Schlag PTI that proves © under the second case of the decay assumption 
(A). 
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Note that, interpolating with the L 2 -bound \\e~ itH °F ac yA|| L 2 < || yr\\ L z, we have 
Corollary 7 Under the decay assumption (A) on the potential V, we have 



where 1 < p < 2 and - + -7 = 1- 

— 1 — p p 

Furthermore, through aT*T argument, © leads to the class of Strichartz estimates 
(see e.g. (Ml): 

Theorem 8 Under the decay assumption (A) on the potential V, we have 



HSCOPflcVlIz^) < C\\ yr\\ L 2, far all - + - = -, 2 < <? < <*>. (7) 

Cj V z. 



1.4 Structure of the paper 

The rest of the paper is organized as follows. In Section [2l we will provide the 
heuristic of the proof of the Theorem [2] In this aim, we will first announce a new 
Theorem [9] providing the same result as the Theorem [2] but under some more re- 
strictive assumptions. After discussing heuristically the proof of this new theorem, 
we will give the elements to extend the proof to that of the Theorem |2] 

The Theorem[9l will be proved in Section[3j Through a rather simple change in the 
feedback law, we will be able to extend this proof to that of the Theorem [2l This 
will be addressed in Section HI 

As it can be seen, the assumptions of the Theorem [2] may still seem too restrictive. 
However, through some arguments based on the analytic perturbation of linear op- 
erators and the quantum adiabatic theory, we are able to relax significantly these 
assumptions. This will be treated in Section [51 



2 Heuristic of the proof 

From now on, we will assume that the system is initially prepared in a purely dis- 
crete state: 



where <^jisc (resp. S &c ) denotes Range(Pdi sc ) (resp. Range(P ac )). The control task is 
to steer the systems state in the eigenspace corresponding to an eigenfunction 0o of 
the free Hamiltonian. Note that this eigenfunction 0o can be any eigenfunction in 




2 N 



^0 — ^O.disc £ <^disc 5 
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the discrete part of the spectrum and does not have to be the ground state. During 
the control process the system might and will cross the continuum <? ac . 

Following the stabilization results for the finite dimensional systems H35I9L a first 
approach for this control problem might be to consider the simple Lyapunov func- 
tion 

y(y) = i-\(y\fo}\ 2 . 

The fact that *P and 0o are both normalized, together with the Cauchy-Schwartz 
inequality, ensures that > 0. The feedback law will be given by ||5): 

u(^) = z((^\h){h\y)), (8) 

where 3 denotes the imaginary part of a complex number. A deep analysis based 
on LaSalle type arguments shows that with such a feedback strategy, one can not 
avoid phenomenons like mass lost at infinity. The population of the state 0o will 
surely keep increasing during the evolution. But, in order to be able to apply the 
LaSalle invariance principle for such infinite dimensional system, one needs to en- 
sure the pre-compactness of the trajectories in L 2 (M. N , C). In the particular case, of 
the Schrodinger equation with the decaying potential, considered in this paper and 
with the feedback law ([8]), one can not even hope to have such a pre-compactness 
result. Indeed, as it has been said before, while the population of the state 0o keeps 
increasing through the application of the feedback law ([8]), during this same period 
some of the population might go through the continuous part of the spectrum. This 
population has then the possibility to disperse rapidly (cf. Subsection 11.31) and so 
we might have some un-controlled part of the L 2 -norm which will be lost at infinity. 

The approach of this paper consists in avoiding the population to go through the 
continuum while stabilizing the state around the target state 0q. So, we consider a 
Lyapunov function which encodes these both tasks: 

M 

:= 1 - (1 - e) £ | (V | (j)j) | 2 - e| | <h) | 2 , (9) 

j=o 

where < £ <C 1 is a small positive constant. Such a Lyapunov function clearly 
verifies: 

and r e (^) =0^1(^100) | = 1. (10) 

Here still, we have used the fact that *P and (j)fs are all normalized in L 2 (M. N ,C). 
Moreover, as the system is initially prepared in the discrete part of the spectrum, 
and as | (¥ | 0o) I > 0, 

%{V ) = 1 - (1 - e) - e| OPo I fo) I < e. (11) 

This Lyapunov function clearly encodes two tasks: 1- it prevents the L 2 -mass lost 
through the dispersion of the absolutely continuous population; 2- it privileges the 
increase of the population in the eigenfunction 0q. 
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By a simple computation we have, 

d M 

= - u (t)[(l - e) £ 3«/iV | 0;> <fc | + £3((m v I' I 0o> <4>o I *»]. 

(12) 

A natural choice is therefore to consider the feedback law: 

M 

u{*¥) = ue(?) := c[(l -e) £ 3(0^ I ^> (fy I ^)) +£3^ | 0o> (fo I *))], 

(13) 

where c > is a positive constant. Such a feedback law clearly ensures the decrease 
of the Lyapunov function y e . Looking at the structure of the Lyapunov function ©, 
the feedback law (TT3T) penalizes strongly exiting from the discrete part of the spec- 
trum. Actually, as i / £ ( x i , o) < £, the decrease in this Lyapunov function ensures that 
the population in the discrete part of the spectrum will always remain more than 
1 — £. Therefore in the worst case, we will only have an £ L 2 -norm which will be 
lost by dispersing in the continuum. 

At the same time, this feedback law (TT3T) slightly encourages the increase in the 
population of the target state 0q. It remains therefore to check whether this increase 
actually provides some kind of convergence toward this eigenfunction or not. This 
will be addressed in Section [3l where we prove the following Theorem: 

Theorem 9 Consider the Schrddinger equation (0Q)- (0. Assume the space dimen- 
sion N > 2. We suppose the potential V(x) to satisfy the decay assumption (A) 
and we take jl e L 2N ~(M. N ) r\L°°(M. N ). We assume the discrete spectrum o c u sc of 
Hq = — A + V(x) to be non-empty. We consider moreover the following assump- 
tions: 

Al = Ylf=o a j$j where {§j}J =Q are different normalized eigenfunctions in the 
discrete spectrum of Ho. 

A2 the coefficient (Xq corresponding to the population of the eigenfunction 0o in 
the initial condition ^q is non-zero: (Xq ^ 0. 

A3 the Hamiltonian Hq admits non-degenerate transitions: Xj l — ^ A /2 — Xk 2 
for (ji,ki) 7^ (72,^2) andwhere {Xj}jL Q are different eigenvalues of the Hamil- 
tonian Hq; 

A4 the interaction Hamiltonian /i(jc) ensures simple transitions between all eigen- 
functions ofHQ.- 

(/%|fe)^0 V/Vfce{0,i,...,M}. 

Then for any £ > 0, applying the feedback law u{t) = u £ {^{t)) given by (TT3T) . the 
closed-loop system admits a unique weak solution in C ([0, r],§)flC 1 ([0, r],//~ 2 (R' 
Moreover the state of the system ends up reaching a population more than (1 — e) 
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in the eigenfunction 0o (approximate stabilization): 

liminf | | fo(x)) | 2 > 1 - £. 

If, moreover, the multiplication by fj.(x) defines a bounded operator over H (Mr), 
then W is a strong solution, i.e. ¥ G C°([0, T],H 2 (R N , C)) DC^fO, r],L 2 (M w , C)). 

Remark 10 77jzs theorem admits some more restrictive assumptions with respect 
to the Theorem ©. In fact, we remove the ID case and we assume the interaction 
Hamiltonian /i to be in a smaller space L 2N ~ (R N ) n L°°(R N ), where 

L 2N -(R N )= |J If, 

2<p<2N 

and therefore, L 2N ~(R N ) C ^(R N ). 

The Theorem [9] will be proved by studying the L 2 -weak limit of *F(f) for t — > °°. 
Namely, let (? n )neN be an increasing sequence of positive real numbers such that 
— > oo. Since || v I / (?„) || L 2 = 1, there exists *Poo G L^R^C) such that, up to a sub- 
sequence, , P(? n ) — ^ weakly in L 2 (M N ,C). Furthermore, through the dispersive 
estimates of the Subsection ll.3[ we will provide a strong convergence result with re- 
spect to the semi-norm = max(||P ( jiscty / 1lL 2 ? HmvIIl 2 )- Through such a strong 
convergence and the Assumptions A3 and A4 of the Theorem, we will prove that 
= /30o> where /3 G C and |j8| < 1. Through some further investigations and ap- 
plying the Assumption Al we will be able to show that |/3| 2 > 1 — £ and this will 
finish the proof of the Theorem |9j 

A deep study of the proof of the Theorem [91 shows that the new restrictions (with 
respect to the Theorem [2]) may be removed if we could ensure the belonging of 
the feedback law w e to the space L l t +8 for 8 G (0, 1]. In fact, the feedback 

law (fT3l only belongs to the space L 2 as ^f- = —\u\. However, we may improve 
this through the following change of the feedback law: 

« £)a (f)= c /(>p)i/(i')r, (U) 

where 

M 

/(¥) := [(1 -e) £3«/i^ I fc) (<{>j I +£3((^ | 0o) (0o I *»], 
and a > and c > 0. This choice of the feedback law implies 

2+q 

and therefore M ea G L/ +a . As a — > °° this ensures that the feedback law w ea be- 
longs to L/ +5 for any 5 G (0, 1]. 
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3 Proof of the Theorem E 



We proceed the proof of the Theorem|9]in 3 steps: 1- we prove the well-posedness 
of the closed-loop system; 2- we prove the existence of an asymptotic regime in 
some appropriate Hilbert space and we characterize the weak co-limit set, i.e. the set 
of the functions 1//00 in L 2 (R N ) such that there exists a sequence of times {t n }™ =l /* 
°° such that ^(tn) — ^ weakly in L 2 ; 3- we finish the proof of the theorem through 
the application of the assumptions Al through A4. 



3.1 Solutions of the Cauchy problem 

Proposition 11 Let £ > and G S. There exists a unique weak solution 
of ©-© with the the feedback law u(t) = MeO^f )) given by (fT3l) . i.e. *P G C°(IR + , S) 
nC l (R + ,H- 2 (R N X)), © holds in H- 2 (R N X)far every t G R + and the equal- 
ity © holds in S. 

If, moreover, G H 2 (R n , C) and multiplication by fj,(x) defines a bounded oper- 
ator over H 2 (R N ), then ¥ is a strong solution, i.e. Y G C°([0, T],H 2 (R N ',€)) n 
C 1 ([0,r],L 2 (R JV ,C)). 

Proo/ Let M G N* be the number of bound states of H = - A + V (x) and T > 
such that 



In order to build solution on [0, T], we apply the Banach fixed-point Theorem to 
the following map 

0: c°([o,r],s) -> c°([o,r],s) 

where ¥ is the solution of ®-0 with u(t) = u E (E, (t)). 

The map © is well defined and maps C°([0,r],§) into itself. Indeed, when § G 
C°([0, r],§), w : f i— > m £ (^(?)) is continuous and thus the Proposition \T\ ensures 
the existence of a unique weak solution V P. Notice that the map takes values 

inc°([o,r],§))nc 1 ([o,r],// (0 2 ). 

Let us prove that is a contraction of C°([0,r],S). Let ^ G C°([0,r],§), uj := 
u E (£j), *¥j := ©(£,-), for j = 1, 2 and A := ^ - W 2 . We have 



2(M+l)c|| J u||i~^ M + 1 )*Hi~ r <l. 



(15) 



A(f ) = i f e 



Jo 



i( - t -^ Ho [u 1 n(x)A(s) + (u 1 -U2)n(x) x ¥ 2 (s))ds. 
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Thanks to (fT3l . we have ||w/||z,°°(o,r) ^ (Af + 1)c||/i||l~ for j = 1,2 and ||vi — 
V2||z,-(o,r) < 2 ( M + - fe|| c o([o,r],L2)- Thus 

< 1)^11^11^1^(^)11^2 + 2(^+1)^11/^11^11^! - <g 2 || c o([o,^],z.2)^- 

(16) 

Therefore, the Gronwall Lemma implies 

l|A(0llco ( [<W2) ^2(M+l)c||M||i~r e ( M+1 )*Hi^ ||^2|| c o([o.r],^ 

and so (fl5T) ensures that is a contraction of the Banach space C ( [0, 7] , §) . There- 
fore, there exists a fixed point \F G C°([0,7],§) such that ©(^ = Since 
takes values in C° ( [0, 7] , §) n C 1 ( [0, 7] , H~ 2 (R N , C) ) , necessarily ¥ belongs to this 
space, thus, it is a weak solution of CO)-© on [0, 7]. 

If, moreover, £ # 2 (K^, C) and multiplication by ju(x) defines a bounded oper- 
ator over H 2 (M. N then applying the Proposition \T\ the map © takes values in 
C° ( [0, 7] , H 2 (R N , C) ) D C 1 ( [0, 7] , L 2 (R n , C) ) thus ¥ belongs to this space and it is 
a strong solution. 

Finally, we have introduced a time 7 > and, for every ^¥q G §, we have built 
a weak solution W G C°([0,r],S) of ©-© on [0,7]. Thus, for a given initial 
condition G S, we can apply this result on [0,7], [7,27], [27,37] etc. This 
proves the existence and uniqueness of a global weak solution for the closed-loop 
system. □ 

Note that, by Assumption Al, the initial state is spanned by the exponentially 
decaying bound states and therefore G Sfl H (M. N ). This, together with the 
Proposition ITTl terminates the proof of the well-posedness part of the Theorem^ 

3.2 Weak co-limit set 

Before studying the weak w-limit set of the closed-loop system, let us announce 
two simple and two rather complicated Lemmas that we will need to characterize 
this asymptotic regime. 

Lemma 12 The feedback law u = u £ (*P) defined by (fT3l) is a member ofL 2 (R + , K). 
In particular, for any J > there exists Ty > large enough such that: 

/*oo 

/ \u E (V(s))\ 2 ds<y. 

J Ty 

Proof. By definition, we have ^ = — ^ |w e ( x P) | 2 . The Lyapunov function 
being a decreasing non-negative function, there exists a positive constant v such 
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that %{^{t)) \ v > 0. Therefore, we have 

\ \u £ (V(t))\ 2 dt = -c ^ = c(r e (^ )-v)<oo. 
Jo Jo at 

□ 

Lemma 13 Let ^(t) denote the weak (or strong) solution of the closed-loop sys- 
tem. There exists a sequence of times (t n )™ =l f oo and some function E L (R , C) 
(with 1 1 *Poo 1 1 ^2 <l) such that: 

*¥„ weakly in L 2 (R N , C), (17) 
^disMtn) -> Prf,A strongly in L 2 (R N , C). (18) 

Proof. The solution \F belonging to C (M + , §), we have ||^2 = 1 and therefore 
the existence of a subsequence Z 1 00 and G L 2 such that (fT7l) holds true is 
trivial. Moreover, 

||^oo||i,2 <liminf||V(r n )|| L 2 = 1. 

n— »o° 

The key to the proof of (fl"8~l) . is in the fact that S^ c = Range(Pdi sc ) is finite dimen- 
sional. Indeed, the weak convergence (fT7l) implies 

| <j)j) -> | , ; = 0,1,---,M, 

and therefore 

M M 

Pdisc^C^) = £ I 0; ) 0; -> £ | (j>j) ^ = Pdisc^oo strongly in L 2 . 

7=0 7=0 

□ 

Lemma 14 Le/; *P(?) denote the weak (or strong) solution of the closed-loop sys- 
tem. Consider a sequence of times (?n)^ = j /* 00 some strictly positive time 
constant x > 0. We have 

S(T)F ac V(t n ) -> ^rangfy m L 2 2 (M^, C) (19) 

Proof. We have the Duhamel's formula: 

S(T)Y(f„) =S(t n + T)V + - f H u e Q¥(s)) S{t n + x-s)n( X yy(s)ds. 

i Jo 

and therefore 

||/iS(T)Pac¥(k)||z2 < ||/x5(? n + T)Pac^0||L2 

+ [ n \u e (^(s))\ \\il(x)S(t n + T-s)F ac ii(x)y(s)\\ L 2ds, (20) 
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where we have applied the fact that the semigroup operator S(t) of Hq commutes 
with the eigenprojection operator P ac of the same Hamiltonian. 

We know by the assumption of the Theorem [9] on /i that /i G Z/(R ) where p G 
[2, 2iV). Applying the Holder inequality, we have 

\\n(x) S(0PacV^lli3 < M L p\\S(t)F ac¥ \\ L ^ + ~ = \ (2D 

where \{feL 2 (R N ). 

Moreover, applying the dispersive estimate of the Corollary Ul we have 

||5(0Pac^lLl < I'f'llvlU ~ + ^ = h (22) 

q q 

for Iff G L q ' nL 2 . 

Let us apply these estimates (|2T|) and (|22~1) to the inequality d20l) . For the first term 
in (|2Q|) . we have 

||/x5(^+T)Pac^0||tf < \tn + A~h\^\\A^t^ ( 23 ) 

X x L x 

where we have used the fact that G L q ' as it is a linear combination of the bound 
states and therefore decaying exponentially ||2). For the second term, we have 

\\nS(t n +T-s)F ac n(xMs)\\ L 2 < \t„ + T-s\ HImII^II^WIL,' 

x *■ L. x 

_N 9 _N - 

<\tn + T-s\ p ||mII^II*(*)IIl2 = \t n + T!-s\ p w\\jp- (24) 



Here, to obtain the second line from the first one, we have applied a holder inequal- 



ity noting that X = \ + ~. 

Furthermore, for any 7 > taking t n > T r (where Ty is given by Lemma O, we 
have 

u e (W(s))\ \\fi(x) S(t n + T-s)F ac Li(x)V(s)\\ L 2d S = 







T r, 







\u e {y(s))\ \\n(x) S(t„ + z-s)F ac n(x)V(s)\\ I 2ds+ 

\u £ (V(s))\ \\n(x)S(t n + z- S )F ac n(x)V(s)\\ L 2ds. (25) 
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Inserting the estimate ((241) in the first integral of (1251) . we have 



\u e {y(s))\ \\n(x) S(t n + T- S )F^(x)^( S )\\ L 2ds < 

/ r T r, ,_av \ 1 2 



< 



\/P o / 2N ~P 2N-p \ 1/2 

7=^ \Wem))\\LjM% [\tn + T-T r \ P -\t n + x\ ~ ) , (26) 

y/\2N — p\ V 7 

where we have applied the Cauchy-Schwartz inequality. Note, in particular that, 
p being strictly less than 2N, ~ p is strictly positive, and therefore the above 
integral (1261 ) tends to as t n — > °°. 

Applying once again the Cauchy-Schwartz inequality, this time for the second in- 
tegral in ([251) , we have 



\u e (y(s))\ \\n(x) S(t n + x - s)F ac ii(x)y(s) \\ L 2ds < 

V2 / rtn 2N N 12 



\u e (¥(t))\*dt J (J (t n + z-s)~~pdsj < 

\/P i/o o / 2N ~P 2W-p \ 1/2 

^^r'^llM^^r— -fe+T-iy-— ) , (27) 

where, we have used the fact that by definition of Ty, \u £ {^{t)) \ 2 dt < y. In 
particular, this implies 



liminf / \u e {y(s))\ \\n(x) S(t n + t - s)F ac ^(x)^(s)\\ L 2d S 



tn^>°° JTy 



t^mhtf 11 ^* 3 *- (28) 



Gathering (|26l) and (|28l) . we have shown 



liminf /" \u E {*¥{s))\ \\n(x) S(t n + T-s)F ac ^{x)f(s)\\ L ids 

tn-* 00 JO x 



-^m^/ ,2Ml;lzr ^- (29) 



Note, however, that we can choose the constant y > as small as we want and 
therefore we have: 

lim f" \u £ (y(s))\ \\n(x) S(t n + T-s)F ac n(x)y(s)\\ L 2ds = 0. (30) 

f«^°°J0 x 
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This, together with (|23l) . finishes the proof of Lemma [141 and we have 



lim ||/xS(T)P ac ^«)|| i? = 0. (31) 

□ 

Applying the above Lemmas, we have the following Lemma, proving the continuity 
of the solution of the closed-loop system with respect to its initial state in the L^ isc - 
topology. 

Lemma 15 Let *P(?) denote the weak (or strong) solution of the closed-loop sys- 
tem. Consider the time sequence {t n }™ =l /* °° and the weak limit as in Lemma\T3\ 
and define W^disc = ^disc^™- Consider the two closed-loop systems 

ij^ n = -A*i> n + V(x)*i> n + u E (y n )n(x)y ni *i> n \ t=0 = V(t n ), (32) 

4* = -AV + V(x)V + u E (V)ii(x)V, n =0 = = V-jisc- (33) 

dt 

We have, for any x > 0, that 

P 'disc^n ( T) -> W disc V( t) strongly in L 2 as n -> oo. (34) 

Proof. In this aim, we consider a stronger semi-norm than L^ isc defined by \\y\\j%> = 
maxdlPdiscV^IlL 2 ' IIMVIIl 2 )- Note however that this semi-norm is weaker than the 
L 2 (M iV )-norm 

HPdiscvy < II vita and ll^vlta < 11/%-llvlta, 

and therefore 

HvH^f < k|IvIIl2(r^ where K = max(l,||/%-). (35) 
It is clear that this is enough to prove 

||¥„(T)-<F(T)||jr->0 asn^°°. (36) 
We have by the Duhamel's formula 

¥„(T)=S(T)¥(fc) + T /'^^Wl^-^W^fv)^ 

«F(t) = 5(t)^ + - [ T u £ (V(s))S(t - s)^i(x)^(s : x)ds : 
I Jo 
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Noting by 5¥ w (t) = ^(t) -^(t), we have 



8y n (T)=S(T)(y(t n )-y ) + - f u e (y n (s))S(T- S )n(x)8y n (s)ds 

i Jo 

+ - f \u E ^ n ( S ))-u e ^(s))]s(t-s)^x) < ¥(s)ds. 
i Jo L -I 



This implies 



\SVn(T)\\je<\\S{T)^{tn)-1fo)\\je 



+ kJ \u e {W n {s))\\\S(T-s)n(x)8V n {s 

'■X 



ds 



+ K 







u e (y n (s))-u e (y(s)) \\S(T-s)n(x)W(s)\\ L 2 m ds, (37) 



where we have applied the inequality (1331) . 

Furthermore, noting that S(t) induces an isometry over the space L (M , C), we 
have 



\\S(T-sMx)8^n(s)\\ L 2 m = \Mx)SV n (s)\\ L 2 {RN) < \\SV n (s)\\jF, 
\\S(t-s)h(x)*¥(s)\\ L 2 {rN) = \\n(x)*¥(s)\\ L 2 {RN) < \\n\\ L ~\\y(s)\\i? < 



where, for the second line, we have also applied 



\m*)\y = \mo)\y <n*-i^<i. 



Inserting the above inequalities in (1371) . we have 



\\8^ n (T)\\^<\\S(T)(^(t n )-^ Q )\\^ + K\\u e \\ LT [ V \\8V n {s)\\jeds 

Jo 



u e (*¥ n (s))-u e (V(s)) 



ds. (38) 



Note, in particular that, by the definition of the feedback law u £ , ||Mg||zr < c(M + 
\u». Let us study the second line of (1381) . We have 



u E (y n (s))-u e (y(s)) 



< 



M 



+ ce 



00 I 
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and for all j E {0,1,..., Af} 



< 

< 

llM^ r »ll^(RV) + ||M||L-||Pdisc5 , i r n(*)||^ CRw ) < (l + llMlMll^n^lU, 

where, for the last inequality, we have applied the Cauchy-Schwartz inequality, the 
facts that 11^0)11/2 = 1 and ||¥(j)|| L 2 < 1, and that | (i/a \ tyj) \ < ||Pdisc V^IL 2 - 

The above inequality, together with (1381) , implies 

\\8W n (x)\\^<\\S{x)^{t n )-W Q )\\^ 

+ cK||/i|| L -(M + 2+||/x|| L -) [ T \\8V n (s)\\^ds. (39) 

Jo 

Applying the Gronwall Lemma to (1391 , one only needs to prove 

\\S(x)(W(t n )-%)\\^^0 asn^oo. (40) 
As a first step, we clearly have 

HPdisc^TX^) -Vo)\\i2 = ||Pdisc*(?«) -* |Il2 - 0, (41) 

where we have used the fact that the semigroup S(x) induces an isometry on L 2 (M. N ) , 
and that the projection operator Pdi sc commutes with the evolution operator S(x). 

Moreover applying the fact that, P disc + P ac = Id L i, we have 

Applying ((4Tt). the first term, ||5"(t) (Fdisc^C^) -^o) IIl 2 converges toward when- 
ever n —> oo. Moreover, applying the Lemma [141 

\\nS(x)F^(t n )\\ L 2^0 asn^oo. 

Therefore, 

\\nS(x){y(t n ) - W ) \\ L 2 - as n - oo. (42) 

The two limits (|4T1) and (1421) imply the limit (1401) and therefore finish the proof of 
the Lemma [T51 □ 

We are now ready to characterize the weak co-limit set. 

Proposition 16 Let *F(f) denote the weak (or strong) solution of the closed-loop 
system. Assume for a sequence (t n )n=i S °° of times that *P(? n ) — ^ & L(M N , C) 
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weakly in L 2 (R W , C) (with H^PooH^ < 1). Define ^oodisc = ^disc^°°- One necessarily 
has 

M 

Uei^dUc) = c[(l -£) £ Sdn^disc I 0;> I *oo )&c )) 

Proof. Consider the Lyapunov function defined in ©. As it is shown in (PT21) . 

the choice (fT3l) of m £ ( v I / ) ensures that the Lyapunov function ^( V I , (?)) is a decreas- 
ing function of time. The Lyapunov function % being a positive function (fTOl) . we 
have 

limr e (^(0) = Tj, (43) 
where 7] > is a positive constant. 

Consider now, the sequence {t n }™ =l f °° of times. The Lyapunov function 
is trivially continuous with respect to i// for the L 2 -weak topology. Therefore, as 
is the weak limit of we have 

W«)= limr e (^„)) = T]. 

Furthermore, noting that the Lyapunov function Y £ only deals with the population 
of the bound states, we have 

r e (v/) = r e (p dlsc v/), 

and therefore 

%{ x V^ c ) = n- (44) 

As in the Lemma [T51 let us consider the closed-loop Schrodinger equation with the 
wavefunction and the initial state ^Vq = ^o^disc- Applying Lemma [T51 for any 
T > 0, we have 

Pdisc^n + t) -> Y(t) Strongly in L 2 (M. N , C) as w ->• °°. 

As the Lyapunov function is continuous with respect to i// for the L^ isc semi- 

norm, we have 

%(V(t n + %))^ %($!(%)) asn^oo. 
But, applying (|43l) . we know that 

r e (»P(f n + T))^77 as/w>o, 

and therefore, 

r e (*(r)) = 7] = r^disc) = r e (*(o)). (45) 
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Thus, the Lyapunov function % remains constant on the closed-loop trajectory of 
*P(f). This, together with CEH) and ([13]), implies 

|-y e (%)) = --^( T )) =0 , (46) 

and therefore by continuity of m £ ( 1 I / (t)) with respect to x and passing to the limit 
at x = 0, we can finish the proof of the Proposition[T6l □ 

3.3 Non-degeneracy assumptions and the proof of Theorem® 

We have now all the elements to finish the proof of the Theorem [9l 

Proposition 17 Let *F(f) denote the weak (or strong) solution of the closed-loop 
system. Consider the sequence (t n )n=i S °°> tne weak limit ^i 1 ^, and its discrete part 
^oo^disc as i n PropositionU6\ Under the assumptions Al through A4 of Theorem® 
we have 

^oo^c = £00, |g| 2 >l-e. (47) 
Proof. Define r\ as in (|43T) . We now, in particular that, 

%{ x ¥ oo ) = r 1 <%(^{0))<e (48) 
where we have applied (flTT) (and therefore the assumptions Al and A2). 
Let us take 

M 

7=0 

Taking the closed-loop system *F(f) as in the proof of the Proposition [T6l we have 
by (|46l) that u £ ( x i , (x)) = 0. Therefore the wavefunction *P(t) evolves freely with 
the Hamiltonian Hq = — A + V(x) and so is given as follows 

~ M 

7=0 

By d46j we have 

~ M 
u(V)=c(l-e) £ g^e'(V^<^|^.) 

+ ce £ g ^ !(Ao ^ )T (ju0j I 0o> = V? > 0. 
7=0 
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The assumption A3 of non-degenerate transitions applies now. As the above rela- 
tion holds true for any T > 0, we can easily see that 

SjGk(nh\ 0j) = O Vj,fce{0,l,..,M}. 

This together with the assumption A4 of simple couplings imply 

gjg k = VMe{0,l,..,M}. 

Thus 

3j e {0,1,...,M} such that q = g 7^ and g = Vfe^j. 

We show that the only possibility for this index j is to be 0. If this is not the case 
(J 7^ 0) taking ¥ = g0 ; - with |g| < 1, 

77 = ^(^Foo) = 1 - (1 -e)|g| 2 > e, 
which is obviously in contradiction with (l48l) . Thus 

¥oo,disC = Sfo 

with I g\ < 1 . Therefore 

%(V 00 ) = l-{l-e)\q\ 2 -e\q\ 2 = l-\q\ 2 . 

Apply once again (|48l ), we have 1 — |g| 2 < £ and so we can finish the proof of the 
Proposition [T71 □ 

Let us now, finish the proof of Theorem [9l 

Proof of Theorem [9| The well-posedness of the closed-loop system has been ad- 
dressed in PropositionQT] In order to prove the approximate stabilization result, let 
us assume that there exists a sequence of times {t n }™ =0 f °° such that 

|Wn)|0o)| 2 <l-e Vn. (49) 

As in Proposition [T71 we can extract from this sequence a subsequence (noted still 
by {t n }™ = o for simplicity sakes) such that 

PdtaW*) L2 ^° ng g> as n -> °o, 
with |g| 2 > 1 — £. This obviously implies 

liminf I I 0o)| 2 > 1-e 

t — >oo 

and is in contradiction with (|49l) . We have therefore finished the proof of the Theo- 
remlU □ 
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4 Proof of the Theorem H 



Let us now get back to the Theorem [2] Comparing to the Theorem |9l the only 
difference is in the fact that, we are also considering the ID case and that the inter- 
action Hamiltonian ji e Jzf HL 00 instead of L 2N ~ PlL 00 for the TheoremUl Therefore 
the only cases remaining to be treated are either the ID case or the cases where 
H e LP(R N )nL°°(R N ) with p > 2N. 

Considering these cases and following the same steps as in the proof of the The- 
orem |9l the only place where we will have a problem to proceed the proof is the 
passage from (T2~6l) and (|27l) to (|28l) and (T2~9l) . Indeed, as IN — p is not strictly posi- 
tive, we can not ensure the convergence towards of the terms in (|26l) and (|27l) . 



A deep study of the estimates (1261) and (|27l) shows that they can be improved if one 
had u £ E L} +s for 5 e (0,1) instead of L 2 as in the proof of the Theorem^ Indeed, 
if one could show that 



u e EL) p for some UJ > 0, (50) 



t 

we could replace the estimates (|26|) and (|27T) with 



o 



\u e (F(s))\ \\n(x) S(t n + x- S )¥ ac ^(x)^( S )\\ L 2ds< 

9 ( f T r \ 1/2 

\\ue(y(t))\\ Li ^- m) Mhy] (tn + V-s) "dsj < 

v/ K0P(r))|| (^^H/ill^ |f„ + T-r y | p -|/„ + t| p 



(51) 



and 



!/2 / rt „ r N \ 1/2 



where 



r p m \ 1 ( f" \ 7 

^ | Me (^(0)|^ _CT Aj (^ + T-j)~djJ < 

\/P 1/9,, ,,9 /, , , ^"A 1//2 

7jfo^«' l l&U T l~-^ + T - r rl~J , (52) 



p-tn(p-N) 
^ N-UJ(p-N)' K } 
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noting that we have applied the Holder inequality and the fact that 

1 1 

+ ■ 



Note that, as p > 2N (this is also true for the ID case as p > 2), G (1,2] and 
therefore there exists some positive QJ > such that — (D > 1 . Furthermore, we 
have 

= p-OJ(p-N) p_ 
^ N-CD(p-N) N' 

and therefore C,N — p is positive. We can thus proceed the proof of the Theorem [2] 
following the same steps as those of the Theorem [9] 



However, it seems that one can not hope to prove an estimate of the form (1501) for 
the feedback law u £ of (TT3T) . We, therefore, need to change the feedback strategy. 
This might be done applying the feedback law (fT4l) . 



The following Proposition clearly implies the Theorem |2] 

Proposition 18 Consider the Schrddinger equation (OQ)- ([2]). We suppose the as- 
sumptions of the Theorem [2] on the potential V (x) and we take n G LP(R N ) n 
L°° (Mr) for some p > 2N. We suppose moreover the assumptions Al through A4 to 
hold true. 

Then for any e > 0, applying the feedback law u(t) = u £ ^ a ( v t '(*)) of CD) with 

p-2N + (D(p-N) „ N 
a = - ; yy ' , 0<G7< 



N-US{p-N) N-p 

the closed-loop system admits a unique weak solution in C ([0, T],S) DC ([0,T], 
H~ 2 (R. N , C)). Moreover the state of the system ends up reaching a population more 
than (1 — e) in the eigenfunction 0o (approximate stabilization): 

liminf | | <^oW)| 2 > 1-e. 

If moreover multiplication by /x(x) defines a bounded operator over H 2 (M. N ), then 
¥ is a strong solution, i.e. ¥ G C°( [0, T] , H 2 (R N , C) ) n C 1 ( [0, T] , L 2 (R^, C) ). 

Proo/ Considering the Lyapunov function % of ©, the choice of the feedback law 
implies 

and therefore proceeding as in the proof of the Lemma [T2l 



2+q g 
T+a r p-JV 



u e , a eLr a =Lr ■ (54) 
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In particular, for 7 > 0, we will chose Ty such that 

/ \u e ( } ¥(s))\p=*~ m ds<Y. (55) 

JTy 

One can then proceed the proof of the Proposition [F71 exactly as in the proof of the 
Theorem|9l replacing only the Lemma[[2]with (|54l ) and (1551 ) and the estimates (|26l ) 
and (EH) by © and {[52]). □ 



5 Relaxations 

As it has been proved in previous sections, the approximate stabilization of a quan- 
tum particle around the bound states of a decaying potential (satisfying the decay 
assumption (A)) may be investigated through explicit feedback laws (fT3l) or (fl4l) . 
The assumptions on the potential V or the interaction Hamiltonian /1 are not so 
restrictive and seem to be satisfied for a large class of physical systems. However, 
the assumptions Al through A4 may seem to be too restrictive. In particular, the 
assumption Al does not allow the approximate stabilization of an initial wavefunc- 
tion with a non-zero population in the absolutely continuous part <c? ac . 

The aim of this Section is to give some ideas to relax these assumptions and to 
consider some more general situations. Some discussions on the assumption Al 
will be addressed in subsection 15.11 Furthermore, a significant relaxation of the 
assumptions A3 and A4 will be addressed in subsection 15.21 Concerning the as- 
sumption A2, we only give the following remark which states that this assumption 
is, actually, not at all restrictive in practice. 

Remark 19 Physically, the assumption A2 in not really restrictive. Indeed, even 
if (^0 I 0o) = 0, a control field in resonance with the natural frequencies of the 
system ( the difference between the eigenvalues corresponding to an eigenfunction 
whose population in the initial state is non-zero and the ground state) will, instan- 
taneously, ensure a non-zero population of the ground state in the wavefunction. 
Then, one can just apply the feedback law of the Theorem\9\or\2\ 

5. 1 Assumption Al 

Before discussing an idea which may result in a significant relaxation of this as- 
sumption, let us provide a remark which states that the result of the Theorem [2] still 
holds true if we relax slightly the assumption Al. 

Remark 20 Consider the Schrddinger equation (OQ)- © with the same assumptions 
on V and ji as in Theorem [2] We consider moreover the assumptions A2 through 
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A4 and we replace Al with 

av \wo\y { w)< T^K^i^o) i 2 . 

The feedback law (fl4l) still ensures the approximate stabilization of the closed-loop, 
one only needs to note that, 

%{V ) < 1 - (1 -e)(l - ^| (¥o | fa) | 2 ) - e| (¥ | fa) | 2 = e, 

where we have applied the assumptions Al ' and A2. The rest of the proof follows 
exactly as in the Theorem® 

Here, we have relaxed the assumption Al by allowing a very small part of the popu- 
lation of the initial state to belong to the continuum. In fact, this allowed continuum 
population is bounded by an 0(e) -proportion of the population in the target state 
fa- 

The assumption Al' of the Remark [20] is still quite restrictive. The question is 
therefore to provide a strategy permitting us to approximately stabilize an impor- 
tant part of the continuum. Note that the controllability of this particular problem 
has never been treated. It seems that one can not in general hope to have a strong 
controllability result. In fact, considering the potentials V and \x of compact sup- 
ports and taking an initial state of support outside supp(V) U supp(/i), it seems that 
an important part of the population may be lost at infinity through the dispersion 
phenomena and this before the controller even has the time to see and to influence 
the state. However, one might be interested to control a part of the continuum. 

Consider for example the potentials V and /i to be negative and of compact sup- 
ports and moreover that supp(/i) C supp(V). Considering the Hamiltonian Hx = 
— A + V + AjUin the strong coupling limit (A — > °°), this Hamiltonian admits more 
and more bound states. One can therefore cover a higher and higher dimensional 
subspace of L (Mr) through the discrete eigenspace of Hx- Assume an initial state 
which has a large population in the continuum of Hq = — A + V but a small 
population in the continuum of Hx for some A > 0. Applying the strategy of the 
Theorem [2l to the free Hamiltonian Hx and the interaction Hamiltonian /i one may 
hope to reach an £-neighborhood of an arbitrary bound state of Hx . Note, in particu- 
lar that while reaching this bound state the control field u(t) has converged towards 
—A. Letting now the control field u(t) ~ —A varying slowly towards zero and ap- 
plying the quantum adiabatic theory (see e.g. [5]) the state of the system will 
follow closely a bound state of the Hamiltonian Hx (A — ► 0). If the target bound 
state of Hx is chosen to be on the analytic branch corresponding to the evolution of 
the ground state of Hq (see e.g. [30J), as the control tends to 0, we may reach the 
e-neighborhood of the desired target fa. This idea of applying the quantum adia- 
batic theory and the large coupling limit to ensure the control of a population in the 
continuum will be explored in future works. 
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5. 2 Assumptions A3 and A4 



After the above discussions on the assumptions Al and A2, let us study the non- 
degeneracy assumptions A3 and A4. Similar assumptions to A3 and A4 have al- 
ready been considered for the stabilization of finite dimensional quantum sys- 
tems D5I9I . When dealing with finite dimensional systems, the assumptions A3 
and A4 are equivalent to the controllability of the linearized system around vari- 
ous eigenstates of the system. For these finite dimensional systems, it was shown, 
in [1331 through the quantum adiabatic theory, and in [9] through the implicit Lya- 
punov control techniques, that the non-degeneracy assumptions can be relaxed sig- 
nificantly. These relaxations have even been applied to the problem of the approx- 
imate stabilization of the quantum particle in an infinite potential well, being an 
infinite dimensional system (see ifTTlO . In this subsection, we will see that such re- 
laxations may also be considered for our control problem of the quantum particle 
in a decaying potential. 

In this aim, we consider the potential V and the interaction Hamiltonian /i both 
to satisfy the decay assumption (A). Similarly to the previous works, we consider 
the family of the perturbed Hamiltonian H a = — A + V(x) + afl(x) with |<r| <C 1 a 
small real constant. The family {H a } is a self-adjoint holomorphic family of type 
(A) in the sense of Kato (see ||30ll . page 375). Thus, the eigenvalues and the bound 
states of H(j are holomorphic functions of o around zero. 

The absence of zero energy eigenstate for Hq = — A + V as been assumed in the 
decay assumption (A) implies the existence of a strictly positive threshold o* such 
that for a evolving in (— a*, a*) the bound states <j)j of Hq stay bound states of 
H(j and do not join the continuum. This ensures that that, the perturbed eigenvalues 
{^i a ,j}f=o °f H a (with Aoj = Ay) are well-defined and remain less than zero. Note 
that, one might have the appearance of new bound states but the bound states of Hq 
will not disappear while considering perturbations of amplitude |c| < a*. We have 
therefore the following Theorem: 

Theorem 21 Consider the Schrddinger equation ©- © with the decay assump- 
tion (A) on both V and jl. We assume moreover that the space dimension N > 2 
and fl E L 2N ~ fl L°°. We consider the assumptions Al and A2 and we replace A3 
and A4 with: 

(A3-A4)' there exists o e (0, a*) such that the non-degeneracy assumptions A3 
andA4 hold for the eigenvalues and the eigenstates of the perturbed Hamiltonian 
H d = -A + V + dn. 

There exists then a feedback law u^), such that the closed-loop system admits a 
unique weak solution and that 

liminf | (¥(f) | 0o)| 2 > 1-e. 
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Remark 22 Roughly speaking the non-degeneracy assumptions A3 andA4 are al- 
ways satisfied unless some kind of symmetry is admitted in the potential. Formally, 
the assumption (A3-A4)' states that if we can break this symmetry through the ad- 
dition of the interaction Hamiltonian, we are still able to ensure the approximate 
stabilization result. 



Remark 23 Applying the same technique as in Proposition\T8\ one can extend the 
result of the Theorem 1271 to the case of dimensions N > 1 and /i G 5£ flL 00 . 

Proof. By the analyticity of the eigenvalues A CT j and (j) a j with respect to o around 
zero, the assumption (A3-A4)' ensures the existence of a strictly positive constant 
G (0, a*) such that the non-degeneracy assumptions A3 and A4 hold true for the 
perturbed Hamiltonians H a with any o in the interval (0, a«). 

Applying once more the analyticity of the bound states {0 J)Cr }jl o with respect to o 
in (0, a*) implies that one can choose <jM G (0, a*) such that 

Mj, a - 0, || L2 < mm(-, 2(M+1)(2 _ £) + 2£ ) 

Vj = 0, 1, • • • ,M and Vcr G (0, c^). (56) 
Now, applying the inequality (l56l) . the assumption A2 implies: 



I (^O I 4>0,a) | 2 > | OF | 00) I 2 " 2 | <^ | 4>0,a - 00 ) | 

>|(^o|0o)| 2 -2||^a-^|| L 2>O, VaG(0,crtt«), (57) 
where we have applied the Cauchy-Schwartz inequality and the fact that 

e 1 

< 



2(M+l)(2-e) + 2e 2' 



Let us consider the Schrodinger equation (E CT ) characterized by the free Hamilto- 
nian H a = — A + V + c/i and the interaction Hamiltonian /! for some o G (0, ct"") . 
Applying the assumption (A3-A4)', the assumptions A3 and A4 hold true for this 
system. Moreover, the inequality (1571) implies the assumption A2 for this system. 
Finally, applying (|56l ), we have {M a + 1 > M + 1 is the number of the bound states 
of H a ) 

M a M 

Packer) (Vo) =1-1 I (^O I fc,cr) | 2 < 1 " £ I (^0 | *j,o) ? 
.7=0 ;=0 
M 

<l-£(|<^ |^a>| 2 -2||^.a-^|| L 2) 

.7=0 

2(M+l) e |pg |fo)| 2 g/2 2 
< 2(Af+l)(2-e)+2e < T3i72l<*°l^l • (58) 
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This implies the assumption Al' or the system (E CT )(having replaced e by e/2). 
Considering therefore the feedback law 

e Ma 

v* m)) := c(l - -) £ | y - >ff > I *(0> 

Z 7=0 

+ c|(/i¥(0|fo,ff)(fo,cx| ¥(*)), c>0, (59) 

and applying the Theorem HI we ensure the approximate stabilization result: 

liminf | Nf{t) | 0o cj) I 2 > 1 - e/2. (60) 

Note that, the feedback law (|59l) means the application of the feedback u := — a + 
v !,ffOP(0) f° r me mam Schrodinger equation (ID)-©- Finally the limit (l60l) implies 

liminf | (¥(0 | O ) I 2 > liminf (\(W{t) | 0o,a) | 2 - 2||fo >tJ - ^oILO 

? >00 f >00 v 4 ' ' 

> l-e/2-e/2 = l-e, (61) 

where once again we have applied (|56l ). □ 
6 Appendix 

This appendix is devoted to the proof of the Proposition [TJ 

Proof of Proposition^ Let G§J > and u G C°([0,r],R). Let 7i G (0,T) be 
such that 

IMU-MIl»(o,t\) < L (62) 

We prove the existence of ¥ G C°([0, T^L 2 ^, C)) such that © holds by apply- 
ing the Banach fixed point theorem to the map 

0:C°([O,r 1 ],L 2 )-,C°([O,r 1 ],L 2 ) 

where *F is the weak solution of 

=H<?¥-u(t)n(x)Z, W(0,x) =V (x) 
i.e. * G C°([0, 71], L 2 ) and satisfies, for every t G [0, 7\], 

W(t) = e- m ° tx ¥ + i f\- m °^u(sMx)^(s)dsmL 2 (M N ,C). 
Jo 
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Notice that takes values in C 1 ([0,ri],//- 2 (M A? ,C)). 

For e C°([0, 7i],L 2 (R",C)), % := ®(&), ¥ 2 := ©(&) we have 

OFi -Y 2 )(0 = * f e- iH ^u{sMx){^ -&){s)ds 
Jo 

thus 

||(Vi-V 2 )(0||^< 11/%-jf -fcllco^r,]^)- 

The assumption (|62|) guarantees that is a contraction of C°([0, T\],L 2 ), thus, 
has a fixed point W e C°([0, Ti],L 2 ). Since takes values in C l ([0, Ti],//" 2 ), then 
*P belongs to this space. Moreover, this function satisfies ©. 

Finally, we have built weak solutions on [0, 7\] for every ^0, and the time T\ does 
not depend on ^q, thus, this gives solutions on [0, T}. 

Let us prove that this solution is continuous with respect to the the initial condition 
¥0, for the L 2 (R N ,C) -topology. Let ^0,^0 e § and 4> the associated weak 
solutions. We have 

|| OP - *) (?) \\ L 2 ^ \\V - 4» || L2 + WnWt- f \u{s) 1 1| OF - *) (s) \\ L zds, 

Jo 

thus Gronwall Lemma gives 

||(V-*)(0lltf<||^O-*0||tf« IWIt " IH ^^. 
This gives the continuity of the weak solutions with respect to the initial conditions. 

Now, let us assume that 6 H 2 (M. N ,C). Take C to be the bound of the multi- 
plication operator /1 over H 2 : i.e. C is a positive constant such that for every <p 6 
H 2 (R N ,C), \\n(p\\ H 2 ^ C|| <p||//2. We consider, then, T 2 > such that C||w||Li(o,r 2 ) 
< 1 . By applying the fixed point theorem on 

2 :c°([o,r 2 ],// 2 )^c o ([o,r 2 ],// 2 ) 

defined by the same expression as 0, and using the uniqueness of the fixed point of 
0, we get that the weak solution is a strong solution. The continuity with respect 
to the initial condition of the strong solution can also be proved applying the same 
arguments as in above. 



Finally, let us justify that the weak solutions take values in S. For e H 2 , the 
solution belongs to C 1 ( [0, T] , L 2 ) D C° ( [0, T] , H 2 ) thus, the following computations 
are justified 



30 



Thus *P(t) E § for every t E [0, T]. 



For £ S, we get the same conclusion thanks to a density argument and the 
continuity for the C°([0, T],L 2 ) -topology of the weak solutions with respect to the 
initial condition. □ 

Acknowledgments : The author thanks K. Beauchard and J-M. Coron for many 
helpful discussions. 
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